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The one-dimensional Holstein model of spinless fermions interacting with dispersionless phonons is 
studied using a new variant of the density matrix renormalization group. By examining various low- 
energy excitations of finite chains, the metal-insulator phase boundary is determined precisely and 
agrees with the predictions of strong coupling theory in the anti-adiabatic regime and is consistent 
with renormalization group arguments in the adiabatic regime. The Luttinger liquid parameters, 
determined by finite-size scaling, are consistent with a Kosterlitz-Thouless transition. 
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The challenge of understanding superconductivity in 
fuUerenes, bismuth oxides, and the high-Tc cuprates has 
renewed interest in models of interacting electrons and 
phonons 0|. Unlike conventional metals these materials 
are not necessarily in the weak-coupling regime where 
perturbation theory can be used or the strong-coupling 
regime in which a polaronic treatment is possible [Q. 
Neither are they necessarily in the adiabatic regime in 
which characteristic phonon energies are much less than 
characteristic electronic energies. This challenge has led 
to numerical studies of the Holstein (or molecular crys- 
tal) model of electrons interacting with dispersionless 
phonons in infinite dimensions, two dimensions, one di- 
mension and on just two sites (see the references in 
The one-dimensional case is important because of the 
wide range of quasi-one-dimensional materials which un- 
dergo a Peierls or chargc-dcnsity-wave (CDW) instability 
due to the clcctron-phonon interaction. Most theoreti- 
cal treatments assume the adiabatic limit and treat the 
phonons in a mean-field approximation. However, it has 
been argued that in many CDW materials the quantum 
lattice fluctuations are important ||^. 

In this Letter we present a study of the one- 
dimensional Holstein model of spinless fermions at half- 
filling using the density matrix renormalization group 
(DMRG). This model is particularly interesting because 
at a finite fermion-phonon coupling there is a quantum 
phase transition from a Luttinger liquid (metallic) phase 
to an insulating phase with CDW long-range order [|,|. 
This illustrates how quantum fluctuations can destroy 
the Peierls state. The Hamiltonian is 
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where Ci destroys a fermion on site i, destroys a lo- 
cal phonon of frequency oj, t is the hopping integral, g 
is the fermion-phonon coupling and a periodic chain of 
N sites is assumed. The phase transition occurs at a 



critical coupling gc separating metallic (0 < g < ^c) and 
CDW insulating phases {g > gc) Hpl. In the strong 



couphng hmit (g^ >> ut) (0) can be mapped onto the 
anisotropic, antiferromagnetic Heisenberg (XXZ) model 
Q which is exactly soluble. The transition occurs at the 
spin isotropy point, is of the Kosterlitz-Thouless (K-T) 
type, and the Luttinger liquid parameters can be found 
in the metallic phase B . 
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FIG. 1. Zero temperature phase diagram of the 
one-dimensional Holstein model of spinless fermions at 
half-filling. For small fermionic-phonon coupling g the sys- 
tem is a Luttinger liquid with parameters that vary with 
the coupling. For large g the system has an energy gap and 
long-range charge-density-wave (CDW) order. The solid di- 
amonds denote the phase boundary from this DMRG study. 
The systematic errors are smaller than the diamonds. The 
results of previous quantum Monte Carlo studies are denoted 
by squares Q and triangles The dotted line is the phase 
boundary from strong coupling theory and the dashed line 
is defined by tj = Amf (where Amf = 8te~'^'^^^^ is the 
mean-field energy gap) and is the approximate location of 
the phase boundary predicted by a two-cutoff renormaliza- 
tion group scheme 

The phase diagram of (|^) over a wide range of adia- 
bacity parameters (0.05 < t/tu < 20) is shown in Fig. 
|l|. A new variant [§ of the DMRG method [0||,0 is 
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used to determine the energy of low- lying excitations to a 
far greater precision than previous quantum Monte Carlo 
(QMC) studies Finite-size scaling (FSS) of a num- 

ber of energy gaps permits the accurate determination of 
(7c and the Luttinger liquid parameters. 

The study of fermion-boson models such as (|l|) by ex- 
act diagonalization or the DMRG presents a challenge 
because there are an infinite number of phonon quantum 
states on each site. Caron and Moukouri have studied the 
XY spin-Peierls and free acoustic phonons models |^ on 
open chains using a conventional DMRG algorithm. The 
simple truncation of the phonon Hilbert space used in 
these calculations can require an excessively large num- 
ber of states, to the extent where the effort expended 
in representing a single site becomes comparable to that 
expended in representing a block. This becomes impor- 
tant when trying to study periodic systems (which are 
more useful for FSS studies) where an extra site is usually 
added to avoid direct interactions between blocks. Jeckel- 
mann and White devised a scheme that maps bosons onto 
fermions which they applied to the polaron problem (a 
single electron interacting with the phonons) in one and 
two dimensions 1^. A more promising method, which 
dramatically reduces the number of states required to 
represent a site, has been used to examine small (6 site) , 
half-filled Holstein systems using exact diagonalization 
pof . We have developed a somewhat similar DMRG al- 
gorithm which is designed to solve periodic systems with 
a large number of degrees of freedom per site. The de- 
tails of the method will be published elsewhere |^ — here 
we concentrate on the results for (|l|). 

The good quantum numbers used are the total fermion 
number N = J^iLi'^l^i, and, for the neutral case (|- 
filled band; N = N/2), the parity (particle- hole) operator 
P : Ci i-^ (— l)'cj;ai i-^ — The energies calculated are 
the ground state energy Eq = EoiN = N/2, P = 1), the 
charge gap Ach = Eo{N = N/2 ±1) - Eq, and the 1- 
and 2- photon gaps (the two lowest neutral excitations 
0) Ai = Ea{N = N/2,P = -1) - Eg and Aa = 
El {N — N/2, P — 1)~ Eq. a number of accuracy checks 
were performed: The DMRG reproduces exact results in 
the non-interacting and strong coupling limits, and the 
DMRG results agree with QMC results for systems of 
up to TV = 16 sites @| within error bars. The DMRG 
accuracy is determined by the parameter m — the number 
of density matrix eigenstates retained per block. Table 
I lists convergence results for Ach, along with the QMC 
results The DMRG errors, being systematic rather 
than statistical, are two to three orders of magnitude 
smaller than the QMC errors. 

Typical FSS plots of the various energy gaps are shown 
in Fig. II for the metallic {g < gc) and insulating {g > gc) 
phases. In the metallic phase the gaps vanish linearly 
with 1/A^ as — > oo, with Ai lying above Ach for large 
N. In the insulating phase Ai lies below Ach and Ach 
and A2 approach non-zero values as N ^ 00 whilst Ai 



TABLE I. Convergence of the charge gap Ach with the 
DMRG truncation parameter m for various system sizes N 
using parameters t = to and g = 1.5uj. QMC results |^ are 
included for comparison. 



m 


N = 4 


iV = 8 


iV = 16 


A = 32 


26 


0.4110 


0.1971 


0.1021 


0.05504 


36 


0.4110 


0.1971 


0.1004 


0.05244 


48 


0.4110 


0.1971 


0.1002 


0.05148 


66 


0.4110 


0.1971 


0.1002 


0.05117 


78 


0.4110 


0.1971 


0.1001 


0.05103 


94 


0.4110 


0.1971 


0.1001 


0.05099 


QMC 


0.416(4) 


0.200(9) 


0.06(3) 





rapidly tends to zero, the state Eo{N — N/2,P = — 1) 
being asymptotically degenerate with the ground state in 
this phase. 

In the QMC studies the critical point gc was deter- 
mined as the point at which an order parameter or the 
charge gap Ach becomes non-zero. However, in a K-T 
transition these quantities behave as Ach ~ e^^'-^^S"^ , 
and there are nonlinear corrections to FSS which make 
the precise determination of gc very difficult by this 
method. Our method of determining gc is inspired by 
work on the frustrated Heisenberg model where the 
transition point was determined by the crossover of sin- 
glet and triplet gaps. It is known that K-T transi- 
tions have a hidden SU{2) symmetry [Q. We hypoth- 
esise that a.t g — gc, the states Eo{N = N/2 ± 1) and 
Eo{N — N/2,P = —1) form a degenerate triplet in the 
thermodynamic limit. Plots of the difference Ai — Ach 
are included in the inset of Fig. ^ for various N. A 
crossover point 5c (A) is defined as the g value at which 
Ai = Ach. gc{N), listed in Table Jll| for various values of 
t/uj, approaches gc as A^ — > 00 |]l4[] ^^The combined errors 
(DMRG truncation, discretization and fitting in g, and 
extrapolation to A = 00) are estimated to be less than 
five precent. 

TABLE II. Convergence of the crossover point gc{N)/uj, 
determined by Ach ~ Ai, with the system size N for various 
hopping parameters t. 



N 


4 


8 


16 


32 


64 


128 


256 


t = O.luj 


2.0878 


2.0911 


2.0920 










t = LU 


1.087 


1.528 


1.591 


1.608 


1.613 






t = IOlu 








2.220 


2.649 


2.765 


2.788 



The resulting phase boundary is shown in Fig. |l|, along 
with the two QMC calculations ^j^, and the result of 
strong coupling theory Q which becomes exact as i — > 0. 
The DMRG results agree well with the strong coupling 
curve for t/uj < 0.2. For large t the results lie close to 

the curve defined by w = Amf = 8te^^*"/^ . This curve 
was predicted to be the approximate phase boundary for 
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FIG. 2. Finite-size scaling of the different energy gaps in 
the (a) Luttinger liquid and (b) insulating phases. The charge 
gap (Ach) and one- and two- photon gaps (Ai and A2) are 
plotted as functions of the inverse lattice size for t — uj 
and (a) g = 1.5ui and (b) g = 1.75a;. Also shown in case (a) is 
Ach as calculated using QMC Q (solid diamonds with error 
bars). In (a) the dashed lines are straight lines through the 
origin, the slope of which can be used to extract the Luttinger 
liquid exponent Kp (see Fig. ^). In (b) the dashed line is a 
guide for the eye. The differences between the two phases are 
seen in: (i) The relative size of Ach and Ai is opposite for 
the two phases, (ii) In the limit N 00 the gaps extrapolate 
to zero (a non-zero value) in the Luttinger liquid (insulating) 
phase. The inset to (a) shows the difference Ai — Ach as a 
function of the coupling for various A''. The critical coupling 
gc is determined as the value at which this difference vanishes 
in the limit A — * 00. 



t > Lo within a two-cutofF renormalization scheme, where 
Amf is the mean- field energy gap jl^ . A saddle-point ex- 
pansion about the mean-field solution suggests that 
there is a first-order transition for uj Amf- 

We next investigate the nature of the transition and 
the Luttinger liquid parameters in the metallic phase 
(0 !i 5 ^ Sc)- For a Luttinger liquid of spinless fermions, 
Eq scales according to £00 — 0, where 

£00 is the bulk ground state energy density and Up is 
the charge velocity. From conformal field theory |lj] 
the scaling forms for the gaps are: Ach ~ '2k''n ^'^^ 
Ai, A2 ^ '^T^'u-p^p ^ where Kp is the correlation exponent. 
The crossover method of determining is equivalent to 
the assumption that Kn = ^ at g = t^c, i-e., the transition 
is of the K-T type . In Fig. || Up (determined from 
the FSS of Eq) and Kp (the values determined from the 
FSS of both Ach and Ai) are shown as functions oi g/uj 
for the case t = O.lw. The Up values agree very well with 
strong coupling theory. The agreement for Kp is not as 
good, due to the presence of nonlinear correction terms 
to the energy gap scaling forms. 
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FIG. 3. Coupling dependence of the Luttinger liquid pa- 
rameters Up (charge velocity) and Kp (correlation exponent) 
for the case t — 0.1a;. The diamonds and stars are the val- 
ues of Kp calculated from the finite-size scaling of the energy 
gaps Ach and Ai, respectively. The solid triangles are the Kp 
values determined from (BV The solid curves are the results 
of strong coupling theory^! . 

A theory for these nonlinear correction terms 
has been developed for the critical case |12|,p0[. 



namely A 



ch 



N 



1 

4K„ 



and A-) 



27rUp 



K„ 



3A 
logW 



logAf ' ■ 

, where A is a constant and Kp 



1/2. By taking the combination 
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ch 



2iiUp 
N 



iK„ 



K„ 



(2) 



the leading nonlinear correction is cancelled at g — gc, 
the next correction being O ^ (logV)^ ) • ^'^^ t — uj and 
g — gc Kp — 0.52 is obtained if (^ is used to determine 
Kp. In comparison, values of 0.59 and 0.42 are obtained 
from the scaling of Ach and Ai, respectively. It might be 
expected that (0) should give better results for Kp around 
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the critical point than the scahng of Ach or Ai. The 
resulting values, plotted in Fig. ^, are in good agreement 
with strong coupling theory. To check the consistency 
of the transition with a K-T transition the value of g at 
which Kp (calculated using (|^)) equals i is listed in Table 
[II. It can be seen that the transition is consistent with 
a K-T transition throughout the phase diagram. 

TABLE III. Transition point (as determined by the 
crossover of Ach and Ai) and g* , the value of g at which 
Kp = i, (where Kp is calculated from (^)) for various hopping 
parameters t. The agreement between g^ and g* is consistent 
with the transition being of the K-T type. 



t/uj 


0.05 


0.1 


0.5 


1 


5 


10 




2.297(2) 


2.093(2) 


1.63(1) 


1.61(1) 


2.21(3) 


2.79(5) 


9* 1^ 


2.299 


2.102 


1.64 


1.62 


2.27 


2.89 



Finally, we consider the question of phonon softening 
and the mixing of phonon and fermion excitations. Fig. 
^ shows the FSS of the energy gaps for a metallic case 
{g < gc) with large hopping t — 5uj. Whilst Ach is lin- 
ear in Ai and A2 are highly nonlinear. This is 
because the lowest fermionic and bosonic, neutral excita- 
tions have the same quantum numbers, those of Ai and 
A2: The non- interacting fermionic gap only becomes 
less than the bare phonon frequency uj for N w ^ w 60 
and thus Ai and A2 are predominantly 1- and 2- phonon 
excitations for small N (flat in only becoming 1- 
and 2- particle-hole excitations (linear in 1/iV) for large 
N. Note that for these parameter values the phonons are 
softened — the renormalised phonon frequency is around 
half the bare phonon frequency to. It would be interest- 
ing to calculate the 1-phonon Green's function to see if 
the phonons soften completely at the transition. The 2- 
phonon Green's function could be used to study phonon 
anharmonicity. 
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FIG. 4. The energy gaps Ach, Ai and A2 as functions of 
l/N for t = 5w in a metallic case g = 2uj. Dashed lines 
are straight lines through the origin. The gaps Ai and A2 
are not linear in 1/A^ because in the adiabatic regime there 
is strong mixing between fermionic and phonon excitations. 
For systems of less than 20 sites the lowest excitations with 
quantum numbers P = ±1 and A'' = N/2 are predominantly 
1- and 2- phonon excitations. 



In conclusion, we have shown that, using a new vari- 
ant of the DMRG, the phase boundary of the one- 
dimensional Holstein model of spinless fermions can be 
accurately determined. The transition is consistent with 
a K-T transition over a wide range of adiabacity. In the 
antiadiabatic limit the phase boundary and Luttinger liq- 
uid parameters agree well with strong coupling theory. 
In the adiabatic limit the phase boundary lies close to 
a curve predicted by renormalization group arguments. 
Challenges that remain include: 1) finding a method of 
cancelling nonlinear corrections to scaling, and hence ac- 
curately calculating the correlation exponent ifp, in the 
whole of the Luttinger liquid regime; 2) developing a the- 
ory of FSS when the conformally invariant field is coupled 
to a dispersionless field with a gap in order to explain the 
nonlinear scaling in Fig. ^; and 3) a detailed investigation 
of phonon softening and anharmonicity. 
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